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Abstract. Weak-strong uniqueness property in the class of finite energy weak so- 
lutions is established for two different compressible liquid crystal systems by the 
method of relative entropy. To overcome the difficulties caused by the molecular 
direction with inhomogeneous Dirichlet boundary condition, new techniques are in- 



£>. \ troduced to build up the relative entropy inequalities 

in 

S Keywords, liquid crystal, compressible hydrodynamic flow, relative entropy, weak-strong 

^O ■ uniqueness. 

en 

o 

1 Introduction 

In this paper, we consider the equations for the compressible flow of liquid crystals. 
We are first concerned with the following simplified version of Ericksen-Leslie system (for 
physical background see El EESl El EESl US] ) : 

p t + div(pu) = 0, 

(p u) t + div(pu ® u) + VP = /_Au - A divfvd Vd - (-| Vd| 2 + F(Vd)I 3 ) V 

d t + u- Vd = e(Ad-f(d)), 

(1.1) 
where p > is the density, u = (ui, U2, U3) is the velocity, d = (di, 62, ds) is the molecular 
direction field, p, A, 6 are positive constants. P(p) is the pressure-density function and 
here we consider the case 

Pip) = *P\ (LI) 

where a > and 7 > 1 are constants. I3 is the 3x3 unit matrix. The term Vd Vd is 
denned by 

Vd©Vd = (Vd) T Vd, 



"Corresponding author. 



where (Vd) T stands for the transpose of the 3x3 matrix Vd. The vector-valued smooth 
function f(d) and the scalar function F(d) are assumed here to satisfy the relation 

f(d) = V d F(d). 

For example, we can choose F(d) and f (d) as 

F(d) = -L(|d| 2 -1) 2 , f(d) = l(|d| 2 -l)d, 
4a a 

where F(d) is the Ginzburg-Landau penalization and o"o > a constant. 

Let O G R 3 be a bounded smooth domain. In this paper, we will consider the following 
initial-boundary conditions : 



and 
where 



(p,pu,d)\ t=0 = (p (x),m (x),d (x)), xett, (1.2) 

u\ dn = 0, d\ dn = d (x) xedQ, (1.3) 

po e L 7 (ft), p > 0; d G L°°(fi) n H 1 ^); 

I |2 

m G L X (Q), m = if po = 0; ^- G L 1 ^). 

Po 

The global existence of weak solutions (p, u, d) to the initial-boundary problem (II. ip - 
( 11. 31) in three-dimension with 7 > 3/2 was obtained by Wang and Yu in [25] and Liu and 
Qing in [23], independently. In addition, when system ( 11. ip is incompressible, Jiang and 
Tan [12] and Liu and Zhang in [22] proved the global weak existence of solutions to the 
flow of nematic liquid crystals for fluids with non-constant density. Based on the existence 
result, Dai et al. in [1] extended the regularity and uniqueness results of Lin and Liu in [18] 
to the systems of nematic liquid crystals with non-constant fluid density. In this paper, 
we are interested in the uniqueness of the weak solution obtained in [25] and [23] . To our 
best knowledge, so far there are very few results concerning uniqueness of weak solutions 
to the initial-boundary value problem (ll.ll) - (ll.3p . 

Second, we consider the following the hydrodynamic flow equation of nematic com- 
pressible liquid crystals : 

p 4 + div(pu) = 0, 

(p u) t + div(pu ® u) + VP = pAu - A div(Vd Vd - -| Vd|%), (1.4) 

d t + u- Vd = #(Ad+|Vd| 2 d), 

where d G S 2 and the other symbols have the same meanings with those in system (11.11) . 
In this situation, the Ossen-Frank energy configuration functional reduces to the Dirichlet 
energy functional. We refer to the readers to consult the recent papers [4] and [23] for 
the derivation of the system ( 11. 4p . For the system ( 11.41) . we are concerned with the same 
initial conditions ( II. 2p but with the following different boundary conditions : 

dd 

u|dQ = 0, Tj-|<9n = xedQ, (1.5) 

where v is the unit outer normal vector of dQ. 
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In contrast with system (jl.ip . from the mathematical point of view, it is much more 
difficult to deal with the nonlinear term |Vd| 2 d appearing in the third equation of Q1.4J) . 
Even for the incompressible flow, there have been no satisfactory results concerning the 
global existence of weak solutions. Recently, Lin et al. in [20] proved the existence of 
global-in-time weak solutions on a bounded smooth domain in IR 2 . For three dimensional 
case, the problem is still open. We should mention that Li and Wang in [15] has established 
its weak-strong uniqueness principle in three dimension, provided that the existence of its 
weak solution is obtained. 

The compressible flow (jl.4p of liquid crystals is much more complicated and hard to 
study mathematically due to the compressibility. For the one-dimensional case, the global 
existence of smooth and weak solutions to the compressible flow of liquid crystals was 
obtained by Ding et al. in j2[ [3]. As for three dimensional case, Huang et al. in [11] and 
Liu and Zhang in [22] established the local existence of a unique strong solution provided 
that the initial data are sufficiently regular and satisfy a natural compatibility condition. 
However, the global existence of weak solution to the compressible flow of liquid crystals in 
multi-dimension is still open. In this paper, we shall establish the weak-strong uniqueness 
property for (11.41) . 

Recently, Feireisl, Jin and Novotny [7J established the weak-strong uniqueness for the 
compressible Navier-Stokes equations. They showed that a weak solution coincides with 
the strong solution with the same initial data, if the strong solution exists. We also refer 
to the readers to consult the recent papers P El El ED] for more weak-strong uniqueness 
results. Motivated by Feireisl et al. [7], we aim to establish the weak-strong uniqueness 
property for two simplified Erichsen-Leslie system (II. ip and (ll.4p . respectively. Our me- 
thod is essentially based on the relative entropy, the modified relative entropy inequality, 
and a Gronwall-type argument. Compared with the existence result of weak solution in 
[231 [23]) where both of them require 7 > 3/2, we are going to make use of the techniques, 
established in [26] . to estimate the remainder 1Z (or 1Z\ (for the definition see (13.101) (or 
(I4.10p )). so as to establish the weak-strong uniqueness property for an improved lower 
bound for any adiabatic exponent 7 > 1. When dealing with the compressible nematic 
liquid crystal flow (II. ip or (II. 4p . the main difficulty lies in the coupling and interaction 
between the velocity filed u and the direction filed d. In particular, we should emphasize 
here that, as far as the weak-strong uniqueness property for the initial-boundary pro- 
blem (ll.ll) - (ll.3p is concerned, more efforts are required to build up the relative entropy 
inequality so as to overcome the effects arising from inhomogeneous boundary condition 
d|an = do(x), which is quite different from the case of the isentropic compressible Navier- 
stokes investigated by Feireisl et al. in [7j. According to the definition of weak solutions 
for the compressible flow of liquid crystals, it is easy to see that 

dGl 2 (0,T;iJ 2 (fl)), 

so the weak solution d actually solves the third equation of (II. ip (or (II. 4p ) in the strong 
sense. Consequently, we can combine the arguments of Feireisl et al. in [7J with the energy 
estimates for parabolic equation to obtain the desired relative entropy inequality. 

The sizes of the positive constants fi, A, and 9 do not play important roles in our 
proofs, we shall therefore assume, for simplicity, that 

n = x = e = i, (1.6) 

throughout this paper. In addition, to simplify the notations, we always set 

A-(VB)C = ((C-V)B)-A, 
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in which A, B, and C are vectors in M 3 . 

The paper of the rest is organized as follows. In the next section, we recall the de- 
finitions of weak and strong solutions to the compressible flow of liquid crystals for two 
kinds of models and state the main results. Section 3 is devoted to the derivation of the 
relative entropy inequality and the proof of Theorem 12.11 Finally, we prove Theorem 12.21 
in Section 4. 



2 Main results 

{p, u, d} is a finite energy weak solution (see [2TJ [21]) to the initial-boundary value 
problem (1TT]1 - (10]1 . if, for any T > 0, 

•p>0, pGL°°(0,T;i^(ft)), uEL 2 (0,T;H^n)), 

d G L°°((0, T) x ft) n L°°(0, T; H\n)) n L 2 (0, T; # 2 (ft)), 

with (p, pu, d)(0, x) = (po(x), m (x), d (x)) for x G ft ; 

• The first equation in (II. ip is replaced by a family of integral identities 

/ p(r, -)(p(T, ■) dx - / p <p(0, -)dx = / / (pd t <p + pu ■ V<p) dxdt (2.1) 

Jn Jn Jo Jn 

for any <p G C l ([0, T] x ft), and any r G [0, T] ; 

• Momentum equations (jl.lft o are satisfied in the sense of distributions, specifically, 

/ pu(r, -)<p(t, •) dx - / p u • <p(0, •) dx (2.2) 

Jn Jn 

= (pu- d t (fi + pu Cg) u : V<p + P(p)divp — Vu : V<p ) dxdt 

+ / / (vdQVd- (i|Vd| 2 + F(d))l 3 ) : V<pdxdt 

for any <p G C l ([0,T] x ft), y?| an = 0, and any r G [0,T]. 

• Equations (ll.ip ^ are replaced by a family of integral identities 



d(r, •) • <p(r, •) dx - / d -<p(0,-)dx (2.3) 

/ / (d ■ dt<p - Vd : V<p - <p • (Vd)u - <p ■ f (d)) dxdt 



for any <p G C^tO.T] x ft), ^| 9n = 0, and any r G [0,T]. 
The energy inequality 



E(t)+ f f (|Vu| 2 + |Ad-f(d)| 2 )dxdt<£(0) (2.4) 

holds for a.e. r G [0,T], where 

E(t) = J (lp|u| 2 + ^p> + i|Vd| 2 + F(d)) dx, 



and 

E(0)= [ (^ + ^ T pU^|Vd | 2 + F(d ))dx. 
Jn ^ 2 Po 1 ~ 1 2 / 

The existence of global-in-time finite energy weak solutions to the initial-boundary 
problem fll . ip ~ fll .3D with the adiabatic exponent 7 > § was established in [231 125], provided 
there exists a positive constant Cq such that d • f(d) > for all |d| > Cq > 0. 

{p, u, d} is called a classical (strong) solution to the initial-boundary problem ( II. ip - 
( TOj) in (0, T) x Q if 

peC^Tjxii), p(t,a;)>p>0 for all (£, x) G (0, T) x fi, 
u, d t u, V 2 u g C([0,T] x tt), d, d t d, Ad e C([0,T] x fi) 

and p, u, d satisfy equation (II. ip . together with the boundary conditions (ll.3p . Observe 
that hypothesis ( 12. 5 p requires the following regularity properties of the initial data : 

p(o,-) = Po ec 1 (n), p >p>o, 

„ ~ _ 2.6 

ii(o, -) = uo6 C 2 (Cl), d(0, •) = do e C 2 (n). 

We are now ready to state the first result of this paper. 

Theorem 2.1. Let Q G M 3 be a bounded domain with a boundary of class C 2+K , k > 0, 
and 7 > 1. Suppose that {p, u, d} is a finite energy weak solution to the initial-boundary 
problem ( fi.il )-( TOl) in (0,T) x Q in the sense specified above, and suppose that {p, u, d} 
is a strong solution emanating from the same initial data $2.6\) . 
Then 

p = p, u = u, d = d. 

In a similar way, we define the finite energy weak solution {pi,Ui,di} to the initial- 
boundary value problem (|1.4j) . (I1.2p . and (II. 5p in the following sense : for any T > 0, 

•Pi>0, Pl eL°°(o,T ; tf(n)), meL^o^H^n)), 

di G L°°((0, T) x Q) n L°°(0, T; tf 1 ^)) n L 2 (0, T; # 2 (fi)), 

with (pi,piui,d!)(0,a;) = (p (x),m (x),d (x)) for a; G fi ; 

• Similar to l|ZIjl - (|2lj]) . the equations flUD hold in £>'((0,T) x ft) 

• The energy inequality 

E 1 (t)+ f f (|Vui| 2 +|Adi + |Vdi| 2 di| 2 )dxdt<£i(0) (2.7) 

holds for a.e. r G [0, T], where 

£iW = J (\pi\ui\ 2 + ^TTTPl + ^iVdil 2 ) da;, 



and 



/" /l |m n | 2 a ., 1 ,„ , l9 \ 

Si0 = / (-^ + — -p2 + -|Vd | 2 )da;. 
Jn K2 Po 7-1 2 y 



Remark 2.1. In the derivation of energy inequality \2.7\ ), we have used the fact that 
|di| = 1 to get 



(fifed! + ui • Vdi) ■ |Vdi| 2 di = ^|Vdi| 2 (fi t |d!| 2 + Ul • V|di| 



0. 



Remark 2.2. We should remark that the global-in time renormalized finite energy weak 
solutions to the initial-boundary value problem \1.4\) , (QUP and M.5\) is still open. For 
one- dimensional case, Ding et al. in |H Ej/ obtained the global-in-time existence of weak 
solutions. 

Similarly, {j5i,Ui,di} is called a classical (strong) solution to the initial-boundary 
value problem (jOj) . ffT2|) and (OJ) in (0,T) x Q if 

Pi £^([0,1] X Q), pi{t,x)>p>0 for all (J, x) 6 (0, T) X ft, 
iii, fifeiii, V 2 Ul g C([0,T] x fi), di, Sfedi, Adi e C([0,T] x ft) 

and pi, Ui, di satisfy equation (J1.4J1 . together with the boundary conditions ( II. 5p . Observe 
that hypothesis ( 12. 8 p requires the following regularity properties of the initial data : 

p 1 (Q,-) = Po eC 1 (ti), Po >p>0, 

ui(0, •) = uo e C 2 (fi), di(0, ■) = d e C 2 (Q). 

We now end up this section with another result of this paper. 

Theorem 2.2. Let Q e M 3 be a bounded domain with a boundary of class C 2+K , k > 0, 
and 7 > 1 . Suppose that {pi, ui, di} is a finite energy weak solution to the initial-boundary 
value problem \1.4\) , U.ty) and M.5\) in (0, T) xQ in the sense specified above, and suppose 
that {pi,iii,di} is a strong solution emanating from the same initial data Ii2.9\) . 
Then 

pi =pi, ui = Ui, di = di. 



3 Proof of Theorem 12.1 



Motivated by the concept of relative entropy in 0, we first define relative entropy 
£ — E f [p, u, d] | [p, u, d] ) , with respect to {p, u, d}, as 

e= f (ip|u-u| 2 + n(p)-n'(p)(p-p)-n(p) + l|vd-vd| 2 )dx, (3.1) 

where 

n(p) = -±- P f. (3.2) 

7-1 

In this section, we are going to deduce a relative entropy inequality satisfied by any 
weak solution to the initial-boundary value problem (jl.ip - (11.3p . To this end, consider a 
triple {p, u, d} of smooth functions, p bounded away from zero in [0,T] x Q, u\an = 0, 
and d|^Q = d (x). In addition, u and d solve the third equation of (II. ip . 

Noticing that the boundary conditions for d and d are inhomogeneous, i.e. d|^ = 
d\an = do(x), we should adapt and modify the arguments in [7] to build up the relative 
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entropy inequality. More regularity of d allows us to make use of energy estimates for para- 
bolic equation, which help us overcome the difficulty due to the inhomogeneous boundary 
conditions. Consequently, combining the arguments in [7] with the energy estimates for 
parabolic equations yields the desired relative entropy inequality. 

To begin with, we take u as a test function in the momentum equation (I2.2p to obtain 

/ pu ■ u(t, •) dx — / po u o • u(0, •) dx (3.3) 

Jn Jn 

= ( pu • d t u + pu (g) u : Vu + P(p) divu — Vu : Vu dxdt 

- I I ((Ad - f (d)) ■ (Vd)u) dxdt. 

Second, we can use the scalar quantity ip = -|u| 2 and ip = W(p), respectively, as test 
functions in the continuity equation (j2.ip to get 

/ -p|u| 2 (r, -)dx = / -p |u| 2 (0,-)dx + / / ( pu • d t u + pu • (Vu)u) dxdt (3.4) 
Jsj 2 in 2 Jo Jn ^ ' 

and 

/ P n'(p)(r, ■) dx = / p n(p)(0, ■) dx + / / (pd t W(p) + pu ■ VlT(p)) dxdt. (3.5) 
Jn Jn Jo Jn v ' 

Since (u, d) solves the third equation of ( 11. ip in the strong sense, it is easy to see that 

d t {d - d) + u ■ Vd - u ■ Vd = (Ad - Ad) - (f (d) - f (d)), a.e.. 
Multiply the above equation by (Ad — Ad) and integrate over Q x (0, r), we have 

-|Vd-Vd| 2 (r,-)dx+ / / |Ad - Ad| 2 dxdt (3.6) 

2 Jo Jn 

= f ^|Vd -Vd(0,-)| 2 dx+ f f (Ad- Ad)-(Vd)udxdt 
Jn 2 Jo Jn 

-If (Ad - Ad) • (Vd)udxdt + f f (f(d) - f(d)) • (Ad - Ad)dxdt. 
Jo Jn Jo Jn 

Finally, multiplying the third equation of (11. ip by Ad — f (d) gives 

/ (-|Vd| 2 + F(d))(r,-)dx+ / / |Ad - f(d)| 2 dxdt (3.7) 

Jn 2 Jq Jq 

= f (^|Vd | 2 + F(d ))dx+ / /(Ad-f(d))-(Vd)udxdt. 

Jn 2 Jq Jq 

Summing up the relations (I3.3l) - ()3.6p with the energy inequality (12.41) . we infer that 
J (Ip|u-u| 2 + n(p)-pn'(p) + i|Vd-Vd| 2 )(r,-)dx (3.8) 



/ 

Jn 



Vu-Vu| 2 dxdt+ / / |Ad-Ad| 2 dxdt 
n./<> Jo Jn 
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< I (i po |u -ii(0,-)r + n(p )-pon / (p(0,-)) + i|Vdo-Vd(0,-)r)da; 



+ p(d t u + u ■ Vu) ■ (u - u) dxdt + / Vu : V(u - u) dxdt 

Jo Jn Jo Jn 

-If {pdttt'{p) + pu • Vn'(p)) dxdt - f f P(p)divudxdt 
Jo Jn Jo Jn 

- f f (Ad - f(d)) ■ (Vd)udxdt + f I (Ad - f(d)) • (Vd)udxdt 
Jo Jn Jo Jn 

+ [[ (Ad - Ad) ■ ((Vd)u - (Vd)u) dxdt 
Jo Jn 

+ f f (Ad - Ad) ■ (f (d) - f (d)) dxdt, 
Jo Jn 

where we have used (13.71) . By virtue of the definition (13.21) of n, it is easy to see that 

u'(p)p - n(p) = p{p) 

and 

/ (pdtU'(p) + pVn'(p) • u + P(p)divu) dx = J d t P{p) dx. 

As a consequence, we deduce from the identities 



P(p)(r,-)dx- / P(p)(0,-)dx= / / <9 t P(p)dxdt 
n Jn Jo Jn 

and (13.81) that the desired relative entropy inequality holds : 

£{t)+ J J |Vu- Vu| 2 dxdt+ / / |Ad-Ad| 2 dxdt (3.9) 

Jo Jn Jo Jn 

< 8(0) + / K(p,u,d,p,u,d)dt, 
Jo 

where 

n = n(p,u,d,p,u,d)-.= n d + n c , (3.10) 

-JZ d := / p(u-u) ■ (<9iU+(Vu)u)dx+ / Vu:V(u-u)dx (3.11) 

Jn Jn 

+ [ {(p- p)d t W(p) + Vn'(p) • (pu - pu)) dx - / divu(P(p) - P(p)) dx, 
Jn Jn 



and 



TZ C := - [ /(Ad-f(d))-(Vd)udxdt+ I [ (Ad-f(d))-(Vd)udxdt (3.12) 
Jo Jn Jo Jn 

+ [ [ (Ad - Ad) • ((Vd)u - (Vd)u) dxdt 
Jo Jn 

+ [ [ (Ad - Ad) • (f (d) - f (d)) dxdt. 
Jo Jn 



c. 



In what follows, we shall finish the proof of Theorem 12.11 by applying the relative 
entropy inequality fl3.9[) to {p, u, d}, where {p, u, d} is a classical (smooth) solution of the 
initial-boundary value problem ( ll.ip . f 12 . 6 1) . and ( II. 5J) , such that 

p(0,-) = p , u(0, -) = u , d(0, -) = d . 

Accordingly, the integrals depending on the initial data on the right-hand side of (13. 9p 
vanish, and we apply a Gronwall type argument to deduce the desired result, namely, 



p = p, u = u, 



d. 



Our purpose is to examine all terms in the remainder ( I3.10|) and to show that they can 
be "absorbed" by the left-hand side of ( I3.9p . 

Compared with [7J , we should remark that the main difficulty comes from the coupling 
and interaction between the velocity field u and the direction field d. Moreover, in the 
context of the weak-strong uniqueness, we only assume the adiabatic exponent 7 > 1, but 
not 7 > 3/2 as in [7J. 

Similar to the proof of Theorem 2.1 in [26] (see also [7]), we use ( 13.1 ip to find that 



n„ 



p(u — u) • ((Vu)(u — u)) da; + 



(P " P) 



Au • (u — u) da; 



(3.13) 



- J divu(p(p) - P'{p){p - p) - P{p)) dx 

- ! ?div(Vd Vd - (-|Vd| 2 + F(d))l 3 ) • (u - u) dx 

/ p(u - u) ■ ((Vu)(u - u)) da; 
Jn 

- J divu(p(p) - P'{p){p - p) - P{p)) dx 



+ 



(p-p) 



Aii - divfvd Vd - (-|Vd| 2 + F(d))l 3 



<> P 

Ad-f(d))-(Vd)(u-u)dx 

(Ad-f(d))-(Vd)(ii-u)dx. 



fu — u) dx 



--■Tl d 



After a tedious but straightforward computation, it follows from ( I3.12p that 



TZ r := n r 



;Ad-f(d))-(Vd)(u-u)dx 



(3.14) 



(Ad - Ad) ■ (f(d) - f(d)) dx + / (Ad - Ad) • (Vd - Vd)udx 

f Ad- (Vd-Vd)(G-u)dx+ f f (d) • ( Vd - Vd) (u - u) dx 
Jn Jn 

(f (d) - f (d)) • Vd(u - u) dx. 



Consequently, it follows from the definitions of TZj and 1Z C that 

TZ(p, u, d, p, u, d) = TZ d + TZ C . 
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(3.15) 



In order to prove the weak-strong uniqueness property, we have to estimate the re- 
mainder TZ. As for IZd- Since the procedures are almost same as that in [26] (see also [?]), 
we just follow [26], list the outlines, and skip the details. In the sequel, we are going to 
focus on the estimation of 1Z C . 

From f )2.6p . it is clear to see that 



/ p(u - u) • ((Vu)(u - u)) dx- f divu (p(p) - P'{p){p - p) - P(p)) dx 

< c||Vu|| LO c (n) £:([p,u,d]|[p,u,d]). 



(3.16) 



Here and hereafter C stands for a generic constant, which may change from line to line. 
Let 



g = g 



; AG, d, Vd, Ad) = Au-div(vd©Vd- (i|Vd| 2 + F(d))l 3 ). 



Obviously, we have 






/ -(p-p)g- (u-u)dx 
Ju P 






= -(p-p)g- (u-u)dx+ / -(p- 

J{§<p<2p} P J{0<P<%} P 


- p) g • (u - 


- u) dx 


+ / -(p-p)g- (u-u)dx. 

J{p>2p} P 







(3.17) 



Similar to that in [26] (see also [7]), we make use of Holder's inequality and Sobolev's 
inequality to show that 



/ I -(P-P)g- (u-u)dx 

{§<p<2p} J{0< P< E } I p 



(3.18) 



<C(5)|||||i3 (a) f^[p,u,d]|[p,u,d]j+5||Vu-Vu||i 2(n) 
for any 5 > 0. On the other hand, noticing that 

n(p)-n'(p)(p-p)-n(p)>Cp^ as p>2p>2p 

and 

P — P I_2 

pa 2 < C, as p > 2p > 2p and 7 > 1, 



PP 
we conclude from the definition of relative entropy £ ( [p, u, d] | [p, u, d] ) that 



/ -(P-P)g- (u-u)dx 

J{p>2p} P 



(3.19) 



{P>2p} 

f 

{p>2p} 



P~ P 



PP 
p-p 



1 
P 2 



pa |u — u| dx 



PP 



1_2\ 2~ Ii~ 1 , 

P 2 2 p 2 g • p 2 |u — u| dx 
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<C||g| 



L°°(n) 



p 7 dx 



P 



u — ul dx 



< C||g||ioo(ft)£f [p,u,d]|[p,u,d]J. 

Next, we continue to estimate TZ C . Note that 

||d||L°o((o,T)xQ), ||d||ioo(( 0i T)xn) < C 
and the fact that f is smooth, we have 

(Ad - Ad) • (f (d) - f (d)) dx 



(3.20) 
(3.21) 



< C 



|d-d| 2 dx 



'a 

<S\\Ad-Ad\\l Ha) + C(8)\\d-d l]L , m 



Ad- Ad| 2 dx 

2 

, 2 (n 

2 



< 8 || Ad - Ad\\i 2(n) + C(8) ||Vd - Vd||£ 2 



(O) 



<<5||Ad-Ad||i 2(Q) + C((5)£^p,u,d]|[p,u,d] 

for any 5 > 0, where we have used Sobolev's inequality. It follows from Holder's inequality 

that 



(Ad - Ad) • (Vd - Vd)udx 



< llulUoo^) || Ad - Ad|| L 2 (n) ||Vd - Vd|| L2(Q) 



(3.22) 



<<5||Ad-Ad|| 2 +C(8) \\u\\l 



- l|2 II VH \7HII 2 

3 (ft) II va ~~ va llL 2 (n) 



< 



n 

Ad| 



L°°(Q) 



<5\\Vu-Vu\\ 2 LHn) 
<8\\Vu-Vu\\ 2 L i {n) 



(CI) 

<5||Ad-Ad||i 2(ft) H-C(5)||u|| 2 oc(n) £([p,u,d]|[p,u,d: 
for any 5 > 0. It is clear that 

Ad • (Vd - Vd)(u - u) dx + f f(d) • (Vd - Vd)(u - u) dx 

l|f(d)|U«(n)) ||u - u||i2(n) II Vd - Vd|| L2(n) 

C{8) (||Ad||L- ( o) + ||f(d)||x- ( o)) 2 IlVd - Vd||| 2(n) 
C{8) (||Ad|U=c (n) + ||f(d)|| Loo(C )) 2 £ ([p,u,d]|[p,u,d] 
Finally, (I3.20p and the fact that f is smooth imply that 

(f(d)-f(d)) -(Vd)(u-u)dx 

L 2 (n) 



(3.23) 



(3.24) 



< C ||Vd|| L oo (Q) ||u - u|| L 2 (n) ||d - d| 



< <J||Vu- Vu||2, 2 



(O) 



< <J||Vu- Vu||| 2(n) 



C(8)(\\Vd\\ L - m y \\Vd-Vd\\ 2 L2m 

C(5)(||Vd|| L » (f ,)) 2 ^([p,u,d]|[p,u,d] 



Summing up relations (I3.9[) ( I3T241) . we conclude that the relative entropy inequality 
yields the desired conclusion 

£([p,u,d]|[p,u,d]Vr) < / h{t)s([p,u,d]\\p,u,d]\(t)dt with some h e L\0,T). 

Thus, Theorem 12.11 immediately follows from Gronwall's inequality. 
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4 Proof of Theorem 12.2 



This section is devoted to proving the weak-strong uniqueness property for the initial- 
boundary value problem f ll .41) . (II. 2j) . and (II. 51) . Compared with the problem discussed in 
Section 3, we would like to point out two points : (i) The term |Vdi| 2 di in the third equa- 
tion of (II. 4p has higher nonlinearity than f(d) in the third equation of (11.10 . so that more 
regularity results should be obtained to overcome the difficulty caused by |Vdi| 2 di. (ii) 
Different from d\gn = do (a:), we now have homogeneous Neumann boundary condition for 

dd 

di at hand, namely, — — \qq = 0, which implies that a modified relative entropy is required 

ov 
to prove the weak-strong uniqueness property for the initial-boundary value problem (11.41) . 

(11. 2p . and (11. 5p . To be precise, we define relative entropy E\ = E\ ( [pi, Ui, di] | [pi, tii, di] j , 
with respect to {p~i, iii, di}, as 

£i = J( ipilui - Gil 2 + (n(p0 - n'(p 1 )(p 1 - p0 - n(p~o) (4.i) 

+ i(|d 1 -d 1 | 2 + |Vd 1 -Vd 1 | 2 ))dx, 

where II is defined by (13. 2p . {p 1; u l5 dx} is the finite energy weak solution to the initial- 
boundary value problem f |1.4|) . (11 .20 . and (11.50 in the sense of section 2, and {pi,tii,di} 

dd 

are smooth functions, pi bounded away from zero in [0, T] x Q, Ui\ga = 0, and — — \qq = 0. 

ov 

In addition, iii and di solve the third equation of ( 11.40 . 

We first establish the relative entropy inequality. For this purpose, take iii as a test 
function in the momentum equation to obtain 

/ /9iUi-Ui(r,-)dx- / p u - ui(0, -)dx (4.2) 

Jn Ju 

— / (piUi • 9ftii + piUi <g) ui : Viii + P(pi)diviii — Vui : Viii J dxdt 

- Adi • (Vdi)tiidxdt 

Jo Jn 

Second, we can use the scalar quantity -|ui| 2 and II' (pi), respectively, as test functions 
in the continuity equation to get 

/ -pi|ui| 2 (r, -)dx = / -p |ui| 2 (0, -)dx+ / / (piiii • 9 t iii H-piGi • (Vui)ui) dxdt 
Jn z Jn z Jo Jn v J 

(4.3) 

and 

I pin , (p 1 )(r, -)dx= [ p n(pi)(o, -)dx+ [ [ ( Pl d t u'{p 1 ) + p lUl ■ vn'(pi)) dxdt. 

Jn Jn Jo Jn v ' 

(4.4) 

Notice that (ui,di) solves the third equation of (11.40 in the strong sense, we have for 
almost everywhere 

$(di - di) + U! ■ Vdi - fli ■ Vdi = (Adi - Adi) + (|Vdi| 2 di - |Vdi| 2 di). 
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Multiplying the above equation by (di— di) and (Adi— Adi), respectively, and integrating 
by parts, we obtain 



we ootam 

[ hd 1 -d 1 \ 2 dx+ [ f \Vd 1 -Vd 1 \ 2 dxdt 
Jn 2 J J n 

f ^|d -d 1 (0,-)| 2 da;+ f / (|Vdi| 2 d! - |Vd! 
Jn * Jo Jn 

pT p 



(4.5) 

2 dx) • (di - di) dxdt 

((Vdi)ui - (Vdi)iii) • (di - di) dxdt 
Jo Jn 

and 

(4.6) 



f -|Vdi- Vdi| 2 dx+ f f | Adi -Adi | 2 dxdt 
Jn 2 Jq Jq 

= f i|Vdo-Vdi(0,.)| 2 da:- f f (|Vdi| 2 di - |Vdi| 2 di) • (Adi - Adi 
Jn * Jo Jn 

+ f f ((Vdi)ui - (Vdi)ui) ■ (Adi - Adi) dxdt 
Jo Jn 

Similarly, we multiply the third equation of (11.4J) by Adi + |Vdi| 2 di to obtain 

/ i|Vdi| 2 (r,-)dx+ f / |Adi + |Vdi| 2 di| 2 
Jn * Jo Jn 

/" 1 r f 



dxdt (4.7) 



-|Vd | 2 dx+ / J Adi- (Vdi)ui dxdt. 



>n l Jo Jn 



< 



Summing up the relations (14 .2ft - (14 .6p with the energy inequality (12.71) . we infer that 

I (i P |ui - ml 2 + n(pi) - pin'(pi) + ifldj - d x | 2 + |Vdi - vdii 2 )) (r, •) dx 

(4.8) 

pT p pT p pT p 

+ / / |Vui - Vdi| 2 dxdt+ / / |Vdi- Vdi| 2 dxdt+ / / | Adi - Ad^ 2 dxdt 
Jo Jn Jo Jn Jo Jn 

J (ipoluo - fii(o, -)l 2 + n(p ) - Pon'(pi(o, ■)) 

+ |(|d - di(0, -)| 2 + |Vd - Vdi(0, -)| 2 )) dx 

+ / / PiiptVii + ui • Viii) • (iii — Ui) dxdt + / / Vui : V(iii — Ui) dxdt 
Jo Jn Jo Jn 

-If (pi^n'(pi) + piui • vn'(pi)) dxdt - / f p(pi)divui dxdt 

Jo Jn Jo Jn 

/ / Adi ■ (Vdi)(ui -Ui) dxdt 
Jo Jn 

f f (Adi - AdO ■ ((Vdi)ui - (Vdi)iii) dxdt 
Jo Jn 

-If (Adi - Adi) ■ (|Vdi| 2 di - |Vdi| 2 di) dxdt 
Jo Jn 

f f (di- di) • (|Vdi| 2 di - |Vdi| 2 di) dxdt 
Jo Jn 



+ 
+ 



JCl 



(di - di) • ((Vdi)ui - (Vdi)iii) dxdt, 



where we have used ( 14. 71) . Thus, we proceed the similar procedures as in the previous 
Section to obtain the desired relative entropy inequality 

£i(r)+ / / |V Ul -Vfli| 2 dxdt+ f I (|Adi- Ad^ + lVdi- Vdi| 2 )da;dt 
Jo isi Jo Jn ^ ' 



< 5i(0) + / 7^i(pi,ui,di,pi,ui,di)dt, 
Jo 



where 



Tli = 7^i(pi,ui,di,pi,ui,di) := TZ ld + TZ lc 



Kid'-= / pi(iii -ui) • ((Vui)(ui -ui))dx 
Jn 

- J divil! (P{ P1 ) - P'(pi)(pi - h) ~ P(Pl)) dx 



(4.9) 
(4.10) 
(4.11) 



+ 



(Pi - Pi) 
<> Pi 



Aui - div( Vdi Vdi - -iVdxl 2 ^ 



(ui — Ui) dx 



and 



TZ lc := / (Adi-(Vdi)-Adi-(Vdi))(ui-ui)da; 
Jn 

+ f (Adi - Adi) • ((Vdi)ui - (Vdi)iii) da; 

;Ad x - Adi) • (|Vdi| 2 di - |Vdi| 2 di) dx 

+ f (dx- di) • (|Vdi| 2 di - |Vdi| 2 di) dx 
Jn 

■ / (di - di) • ((Vdi)ui - (Vdi)iii) dx. 



(4.12) 



Similar to the proof of Theorem 12.11 we shall finish the proof of Theorem 12.21 by ap- 
plying the relative entropy inequality (14. 9 j) to {pi, Ui, di}, where {pi, u 1} di} is a classical 
(smooth) solution of the initial-boundary value problem (jl.4p . (j2.6p . and (II. 5p . such that 

Pi(0,-) = Po, iii(0,-) = Uo» di(0,-) = d . 



As a result, the integrals depending on the initial data on the right-hand side of (14.91) 
vanish, and we apply a Gronwall type argument to deduce the desired result, namely, 

Pi = Px, ui=ui, di = di. 



To complete the proof, we have to examine all terms in the remainder (I4.10p and to show 
that they can be "absorbed" by the left-hand side of (14.91) . 
From ( 12.61) . it is clear to see that 

/ pi(tii - m) ■ ((Viii)(ui - ui)) dx- f divfix (p( Pl ) - P'(pi)(pi - h) ~ P(P±)) &x 
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< C||V(ii||i,cx)(fj)f 1 . 
Let 



(4.13) 



gi = gi(Aui, Vdi, Adi) = ACn - div(Vdi © Vdi - -|Vdi|% 
Obviously, we have 

/ — (pi-pi)gi-(ui-ui)da; (4.14) 

Jn Pi 

= / — (pi - px) gi • (ui - ui) dx + / — (pi -pi)gi • (ui -ui)dx. 

</{0<pi<2pi} Pi J{pl>2pi} Pi 

Similar to (|3.18p . it follows that 



— (Pi - Pi) gi • (ui - ui) dx 

{0<pi<2pi} Pi 



51 112 



<c{8)\m\Uw£i 

Pi 



/ — (Pi-Pi)gi- (iii-ui)da; 

J{ P l>2pi} Pi 



Now we are going to estimate TZi c . From (14.121) . a direct calculation gives 



n lc = n a lc + n 



lc 



where 



K 



a lc = /" (Adi - Adi) • (Vdi - Vdi)iii da; 

Jn 

+ / Adi- (Vdi- Vdi)(iii-ui)dx 



(|Vdi| + |Vdi|)di • (Adi - Adi)(|Vdi| - |Vdi|) dx 



+ / (di-di)-(Vdi)(ui-iii)dx 
'a 



and 






(di - di) • (Vdi - Vdi)tii dx 



+ / (|Vdi| + |Vdi|)di • (di - di)(|Vdi| - |Vdi|) dx 
'a 



|di-di| 2 |Vdi| 2 dx. 



For the first term on the right hand side of (14.181) . we have 



/ (Adi - Adi) • (Vdi - Vdi)ui dx 
Jn 

< ||ui||z,°°(n)||Adi - Adi|| L 2 (n) ||Vdi - Vdi|| L 2 (n) 
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(4.15) 



for any 5 > 0. We proceed in a same way as the derivation of (13.191) to find that 

<C||gi|| LO o ( n)£i. (4.16) 



(4.17) 



(4.18) 



(4.19) 



(4.20) 



C(S)\\u 1 \\ 2 Lx(n) £ 1 , 



< 5 || Ad! - Adi||£ 2(a) + C^lluillioo^HVd! - Vdillia^ 

<6\\Adi-Adi\\%n Q) + C(S)\\i 
for any 5 > 0. In a similar way, we also obtain 

' Adi ■ (Vdi - Vd x )(ui - ui) dx 



< <J||Vui - Vui||^ 



L 2 (n) 



■C(S)\\Ad 1 \\l 00{n) £ 1 , 



(4.21) 



- / (|Vdi| + |Vdi|)di ■ (Adi - Adi)(|Vdi| - |Vdi|) dx 



and 



< 5 || Adi - Adi||^ 2(n) + C(5) (HVdxHioo^ + HVdiH^^JIIdi 



(di — di) • (Vdi)(ui — ui) dx 



I 2 F 

lL°°(n) c i) 



< 5 ||Vui - Vui||i 2(n) + C{8) || Vdi|||„ (n) 5i, 

for any 5 > 0. 

Finally, Cauchy inequality and the definition of the relative entropy E\ give 

(di-di)-(Vdi-Vdi)uidx 

< ||uilU<»(fi)||di - di|| L 2( Q) ||Vdi - Vdi|| L 2 (Q) 

i i 

<l I ~ 1 1 C 2 C 2 

— \\ u l\\L°°(fl) &\ Z>\ 

< ||iii||icx>(n) £i- 
Likewise, we get 

(|Vdi| + IVdxDdx ■ (dx - di)(|Vdi| - |Vdi|) dx 



(4.22) 



(4.23) 



(4.24) 



(4.25) 



and 



< ||di||ioo(n)(||Vdi||ioo(n) + ||Vdi||L<»(n))||di — di|| L 2(^)||Vdi — VdiH^n) 

< ||di||ioo(n)(||Vdi||io (n) + ||Vdi||ioo(f2)j S\ 



|di - di| 2 |Vdi| 2 dx < ||Vdi||*oo(n)||di - di|||a (ft) < ||Vdi|||oo (n) £i. 



(4.26) 



Moreover, by applying the quasilinear equations of parabolic type estimates (see [13 
Chapter VI, Section 3) to the third equation of (jl.4p and a density argument, we have 



|Vdi|| L oo (a) <C. 



(4.27) 



Summing up relations fl4.10p - fl4.27p with the definition of finite energy weak solution 
{pi,ui,di}, we conclude from the relative entropy inequality (14.91) that 

£i([pi,ui,di]|[pi,ui,di]j(r) < / h l (t)£ l ^[p 1 ,u 1 ,d 1 ]\[p 1 ,u 1 ,d 1 ]J(t)dt, 

where hi G L 1 (0,T). Thus, Theorem 12.11 immediately follows from Gronwall's inequality. 
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Remark 4.1. For one- dimensional case, Ding et al. in JUG]/ established the existence of 
global classical solutions and the existence of global weak solutions. Based on the existence 
results, one can directly apply Theorem \2.2 to obtain the weak-strong uniqueness principle 
for one- dimensional compressible liquid crystal system. 
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